Observation of anisotropic diffusion of light in compacted granular porous materials 

Erik Alerstam' Ifl and Tomas Svensson''^'[l| 

'Department of Physics, Lund University, P.O. Box 118, 221 00 Lund, Sweden 
^European Laboratory of Non-linear Spectroscopy (LENS), 50019 Sesto Fiorentino (Florence), Italy 

It is known that compaction of granular matter can lead to anisotropic mechanical properties. Recent work 
has confirmed the link to pore space anisotropy, but the relation between compression, mechanical properties 
and material microstructure remains poorly understood and new diagnostic tools are needed. By studying the 
temporal and spatial characteristics of short optical pulses diffusively transmitted through compacted granular 
materials, we show that powder compaction can also give rise to strongly anisotropic diffusion of light. Inves- 
tigating technologically important materials such as microcrystalline cellulose, lactose and calcium phosphate, 
we report increasing optical anisotropy with compaction force and radial diffusion constants being up to 1.7 
times the longitudinal. This open new and attractive routes to material characterization and investigation of 
compression-induced structural anisotropy. In addition, by revealing inadequacy of isotropic diffusion mod- 
els, our observations also have important implications for quantitative spectroscopy of powder compacts (e.g., 
pharmaceutical tablets). 



The physics of granular media in general, and compression 
and deformation in particular, is utterly complex fP-'?!. Ma- 
terial anisotropy due to compression is one phenomena in this 
context that remains poorly understood. It is well known that 
uniaxial compression of certain granular media can induce 
anisotropy in mechanical properties |5, 6|, and scanning elec- 
tron microscopy, pulsed-gradient stimulated-echo NMR and 
spin-echo small-angle neutron scattering have recently con- 
firmed that the pore structure itself indeed can be anisotropic 
IIT HToI . The fundamental understanding of the interplay be- 
tween compression, microstructure and anisotropy is, how- 
ever, still in its cradle Q. 

In this work, we report that compression of granular mat- 
ter also can give rise to anisotropic diffusion of light. Since 
the phenomenon is linked to pore space anisotropy, light scat- 
tering and diffuse spectroscopy may thus turn out to be a 
valuable, and non-destructive, tool for material characteriza- 
tion and fundamental investigations of compression-induced 
anisotropy. As the materials in which we observe anisotropic 
diffusion of light are of major technological importance, our 
findings also reveal an urgent need for better understanding of 
the microstructure and optics of compacted granular media. In 
particular, our results have important implications for quanti- 
tative spectroscopy of, e.g., compacted powders and pharma- 
ceutical tablets. Anisotropy can, for example, be a compli- 
cation for spectroscopic methods that aim at separating the 
effects of scattering and absorption, methods which typically 
are based on isotropic diffusion models [11 J. 

Anisotropic diffusion of light has previous been observed 
and studied in materials with rather evident anisotropic mi- 
crostructure, including etched porous semiconductors [12J, 
nematic liquid crystals in3tfT6l . fibrous and stretched plas- 
tics |TT, 18], wood fT9l and biological tissues such as mus- 
cle, teeth, bone, brain and skin |20-22|. In contrast, we re- 
port on anisotropic diffusion in materials where presence of 
microstructural anisotropy is far from self-evident, still be- 



ing the subject to active research, and where Ught propagation 
conventionally is assumed to be isotropic. In fact, the radial 
symmetry of the diffusion occurring in the compacted granu- 
lar systems discussed here makes it impossible to detect op- 
tical anisotropy by, e.g., employing steady-state methods and 
searching for spatial asymmetries in reflected or transmitted 
intensity profiles. To reveal the anisotropic diffusion of light, 
we have instead relied on spatially resolved measurements of 
light transport dynamics. The main concept behind our ap- 
proach is illustrated in Fig. [T] 
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FIG. 1 : Measurement principle. In spatially resolved photon time- 
of-flight spectroscopy, anisotropic diffusion manifests itself as an 
anomalous spatial dependence of the arrival and shape of diffusively 
transmitted pulses. When radial diffusion is faster than longitudi- 
nal (Dr > Dz, in terms of diffusion constants), as for the cylin- 
drical compacts involved in this study, pulses detected at non-zero 
radial positions arrive earlier and are narrower than what could be 
expected from looking only at on-axis transmission. In fact, the tem- 
poral shape of the pulse transmitted along the axis of incidence is 
independent of the radial diffusion constant. 
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The granular materials investigated are cylindrical com- 
pacts made by uniaxial compression of powders based on mi- 
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crocrystalline cellulose (MCC), lactose or calcium phosphate. 
All samples have a diameter of 13 mm and a total weight of 
500 mg. The compaction force applied varied between 5 and 
55 kN, and thickness were between 1.9 and 3.5 mm. The 
MCC-samples represent realistic pharmaceutical wet granu- 
lated compacts, and were manufactured from three different 
granule size fractions A, B and C (granules being < 150 [a.m, 
150 — 400 |J.m and > 400 |xm, respectively). Lactose and cal- 
cium phosphate samples were, on the other hand, made from 
pure powders (5-10 [a.m mean particle size). For reference, we 
also make experiments on homogenous isotropic turbid mate- 
rials with scattering properties similar to those of the granu- 
lar samples. These are Spectralon (a commercially available 
porous fluoropolymer), a macroporous sintered alumina ce- 
ramic [23, 24 1, and a Ti02 -based epoxy phantom |25|. 

Light transport is investigated by conducting spatially re- 
solved photon time-of-flight spectroscopy (PTOFS). Short 
(picosecond) pulses of 760 nm light are injected into the 
tablets at the radial position Ts using a 600 /im-core graded- 
index optical fiber, and light transmitted through the sample 
is detected at the radial positions using a second identical 
fiber (radial positions defined relative to the sample center). 
A key parameter in measurements is the radial (transverse) 
source-detector separation, Sr = — rd- The collected light 
consist of pulses that are temporally broadened due to multi- 
ple scattering, and the detection fiber direct these pulses to a 
fast photomultiplier that combined with time-correlated single 
photon counting allows us to resolve them in time. Details on 
the spectroscopic instrument is available in Ref. [25], and the 
experimental configuration is elaborated in Figure |2] 




Detector fiber 'Tj 



FIG. 2: Experimental configuration. The detector fibers runs 
through the center of a cylindrical block. The sample is put on top, 
surrounded by a ring of having an inner diameter matching the sam- 
ple. To reduce stray light, parts are made in black Delrin plastic. 
Part (a) shows arrangement spatially resolved PTOFS. Centering of 
sample and ring sets ~ 0. Sideway translation of the source fiber 
allows tuning of rs . Part (b) shows how transverse translation of sam- 
ple and ring allow simultaneous tuning of Vd and rs, a configurations 
used for investigation of sample homogeneity (keeping 5r = 0). 

To avoid being very close to the edge of the samples (sam- 
ple radii being 6.5 mm), measurements were conducted up 
to Sr — 5 mm. For the most strongly scattering samples, 
however, there was an upper limit to Sr above which the sig- 
nal is so weak so that stray light distorts obtained time-of- 
flight (TOF) histograms. In calcium phosphate samples and 
the porous alumina, for example, light completely change its 
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FIG. 3: The failure of isotropic modeling. For homogenous 
isotropic materials, isotropic diffusion modeling should result in the 
same fitted diffusion constant regardless of the source-detector sep- 
aration used in the measurement. While this is the case for our 
isotropic reference materials (solid dots), it is evident that the com- 
pacted granular samples (open circles) cannot be modeled in the 
same manner For clarity, the graph shows data only from a few of 
the samples (compression force in kN is given within parantheses). 



direction on average each every 5 /im (i.e., the transport mean 
free path of light is around 5 /im). The resulting high reflectiv- 
ity and strong attenuation limits us to measurements up to 8r 
of 3 mm. At larger 8r, measurements are susceptible of stray 
light, e.g. light that exited the sample at lower br and ideally 
should not reach the detector fiber. Analysis of the obtained 
photon TOF histograms are made by employing standard dif- 
fusion models generalized for anisotropic diffusion lfT4l [TSl , 
and using the extrapolated boundary conditions as described 
by Contini et al. |26|. 

The inadequacy of homogenous isotropic diffusion mod- 
els is elucidated in Fig. [3] There, measurements at different 
source and detector separation are evaluated individually us- 
ing the isotropic diffusion model. If diffusion is isotropic, fit- 
ted diffusion constants should be constant (independent of (5r). 
Data from the reference materials are thus in excellent agree- 
ment with isotropic diffusion, and shows the capability of our 
system to deliver consistent data over a wide range of optical 
properties and source detector separations. At the same time, 
it is evident that the light transport in the compacted granular 
samples cannot be explained by homogenous isotropic diffu- 
sion. Note that also for the calcium phosphate data, although 
not easily seen in the figure, fitted diffusion constants are non- 
constant and systematically increasing with 8r (growing from 
0.29 at 8r = to 0.33 mm^/ns at br = 3). 

One possible cause of the mismatch with the model of ho- 
mogenous isotropic diffusion could, of course, be that our 
samples exhibit a spatially varying diffusion constant due to 
macroscopic density heterogeneity. However, careful investi- 
gation of sample homogeneity rules out this possibility. First, 
we investigated the radial sample homogeneity optically by 
making PTOFS measurements with (5r = at different po- 
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FIG. 4: The spatio-temporal picture. When moving from isotropic 
to anisotropic modeling, fits obtained during evaluation of spatially 
resolved PTOFS (i.e., evaluation of Sr series) goes from being clear 
evidence of model inadequacy to becoming virtually perfect. The 
displayed data is from the MCC C50 sample but is representative for 
all samples. For visual clarity, only a subset of the involved spatial 
recordings are shown (Sr — 0, 1, 2, 3, 4 mm). 



sitions of the samples (r^ = from to 5 mm in steps of 
0.5 mm, cf. Fig. For the compacted granular samples, 
coefficients of variation in fitted D was smaller than 3% (sim- 
ilar variation was seen also for the isotropic reference materi- 
als), and dependence on radial position was insignificant. This 
shows that our compacted samples are radially homogenous. 
To verify longitudinal homogeneity, our structures were also 
investigated by X-ray microtomography (/iCT). Since radial 
homogeneity is confirmed by our optical experiments, /iCT 
scans were optimized for longitudinal density mapping. Still, 
there were no signs of longitudinal heterogeneity. 

Convincing evidence that anisotropy is the proper explana- 
tion of why homogenous isotropic diffusion modeling fails is 
presented in Fig. |4] There, simultaneous evaluation of the 
TOF histograms of a 6r series clearly shows that isotropic 
homogenous diffusion is far from being capable of explain- 
ing experimental data, while the anisotropic diffusion model, 
with its two diffusion constants, results in excellent fits. Fig. |4] 
gives only one example, but similar results were observed for 
all compressed granular samples. Even for the calcium phos- 
phate samples, which looking at Fig. [3] may appear close to 
isotropic, a dramatic reduction of residuals are obtained when 
moving to anisotropic diffusion modeling. In comparison, and 
as expected, anisotropic modeling of data acquired from the 
isotropic reference materials result in insignificant improve- 
ment of fits. 

It is important to realize that estimation of Dr requires mea- 
surements at Sr 7^ 0, while can be estimated using only a 
single measurement at Sr = 0. This can be inferred directly 
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FIG. 5: The necessity of spatial information. Typical evolution 
of fitted radial and longitudinal diffusion coefficients (Dr and D^) 
as more and more spatial information is included (i.e., evaluation 
involving all datasets up to an upper Sr). In order to reach good esti- 
mates of both longitudinal and radial diffusion constants, evaluation 
must clearly include measurements made at Sr being on the order 
of the sample thickness. Here, datapoints from three measurement 
repetitions from the MCC A (20kN) sample are shown. 



from anisotropic diffusion theory (cf. Eq. [TJ, and the issue 
is elaborated further in Fig. |5]by showing how estimates of 
diffusion constants evolve as the amount of spatial informa- 
tion increase. While estimation of is very robust, proper 
evaluation of requires inclusion of measurements made 
at (Jj'-values comparable to sample thickness. On the other 
hand, Dr converge to a constant value rather quickly. This 
means that inclusion of very large (5r-values is not necessary 
in general, making it possible to avoid low transmission mea- 
surements where stray light may distort data. Measurements 
made very close to the radial boundary may also be influenced 
by edge effects in material density, and eventually also by the 
breakdown of the infinite slab diffusion-model. In fact, seen 
in Fig. |5j we generally observe a minor change in Dr when 
including measurements at high Sr (when Sr — 5 mm, the 
outer part of the fiber core is only 1.2 mm away from the sam- 
ple boundary). In the following, our evaluation of anisotropy 
is therefore based on simultaneous evaluation of TOF his- 
tograms recorded at Sr not larger than 4 mm (as mentioned 
earlier, the strongly scattering calcium phosphate samples are 
exceptions, and evaluation is limited to Sr not larger than 3 
mm). 

Fig. [6] presents observed anisotropy and longitudinal dif- 
fusion constants for all the materials. Measurement repeti- 
tions (also shown) show that the variation in derived D^, Dr 
and 7 generally is very small. The coefficients of variation 
were between 0.3 and 3.3%, except for the 50 kN calcium 
phosphate which gave a 6.6% percent variation in Dr and 7). 
The anisotropy of the MCC and lactose samples increase with 
pressure. Interestingly, this trend is not shared by the calcium 
phosphate samples. In this context, it is important to note also 
differences in scattering properties. It is clear that the scatter- 
ing strength is largely the same for the calcium phosphates. 
In contrast, scattering decrease with pressure for lactose and 
MCC samples (note, however, the plateau reached at 30 kN 
for MCC samples, after which scattering remains constant). 
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FIG. 6: Anisotropy summary. Longitudinal diffusion constants and 
anisotropy factor 7 = Dr /D^, versus compression force for all sam- 
ples. The general trend is that both anisotropy and diffusion constants 
increase with pressure (the increase in D corresponds to a decrease of 
the scattering strength). MCC samples are shown as solid circles (•), 
plus signs (+) and open circles (o), for the granule size classes A, B 
and C respectively. Lactose samples are given by triangles (A), and 
calcium phosphate samples by the times symbol ( x ). Note that three 
measurement repetitions were done on samples made at the lowest 
and highest compression force, as well as for the intermediate 20 kN 
MCC samples. 



DISCUSSION 

The fact that Hght propagation in compacted granular me- 
dia can be strongly anisotropic represents a major revision of 
our understanding of the optics of such systems. It is there- 
fore important to realize that our findings correlate well with 
recent reports on anisotropic pore structure. For example, 
pore space anisotropy of MCC compacts have been confirmed 
and shown to increase with compression force f8^ TU]. In 
addition, the pore space anisotropy reported in these studies 
occur on length scales that should influence light scattering. 
The fact that radial diffusion is faster than longitudinal is pre- 
sumably related to formation of crack-like pores running per- 
pendicular to the compression axis [7j. On the other hand, 
it should be noted that NMR porosimetry has been applied 



to e.g. lactose and calcium phosphate samples without clear 
signs of anisotropy. However, the major difference in the way 
the the the pore space is probed (and differences in samples) 
makes deeper comparison meaningless. We note, instead, that 
anisotropy of mechanical properties have been observed for 
compacts based on all material investigated here i5]|6l, ren- 
dering our observations highly plausible. In fact, since light 
scattering and its wavelength dependence is extremely sensi- 
tive to microstructure, we believe that spectroscopic optical 
methods may turn out to be very useful and sensitive probes 
of pore space anisotropy. Conventional porosimetry methods 
like mercury intrusion and gas adsorption does not give infor- 
mation of pore dkectionality, and currently employed meth- 
ods that do are subject to important limitations [9 1. 

Besides opening a promising new route for characteriza- 
tion of material anisotropy, our results also have obvious im- 
plications for quantitative spectroscopy of compacts in gen- 
eral. The use of isotropic diffusion models must be recon- 
sidered. Moreover, considering a combination of anisotropic 
effects with the density variations characteristic of more com- 
plex shapes of e.g. pharmaceutical tablets fTJ], it is clear that 
model-based quantitative spectroscopy of powder compacts 
indeed is challenging. The picture gets even more compli- 
cated when considering also the more fundamental issue of 
how light sample complex porous structures Il23ll28l . 

To conclude, our observation of strong anisotropic diffu- 
sion of light in compacted granular media have broad impli- 
cations. As discussed above, the phenomenon can be used 
for material characterization and fundamental investigations 
of compression physics, and also has important implications 
for quantitative spectroscopy of compacts. At the same time, 
reaching full understanding of the relation between compres- 
sion, microstructure and light scattering is challenging, and 
this endeavor will require extensive and interdisciplinary ef- 
forts in the future. 



METHODS 

Compacted granular samples 

The compacted granular samples are based on either mi- 
crocrystalline cellulose (MCC), lactose or calcium phosphate. 
All compacts were manufactured using a standard single- 
punch tablet press with a cylindrical die (all having a 13 mm 
diameter and a total weight of 500 mg). The MCC-based 
samples represent realistic pharmaceutical tablets, contain- 
ing 38.1% MCC (excipient), 31.8% metoprolol (active ingre- 
dient), 12.7% sodium starch glycolate (disintegrant), 11.1% 
lactose (filler), 2.9% colloidal silica, 2.2% povidone (binder) 
1.2% magnesium stearate (lubricant, reduces friction during 
compaction). In order to investigate influence of granule 
sizes, MCC compacts were manufactured from sieved gran- 
ules (granule sizes: (A) < 150 ^m, (B) 150-400 fim, and (C) 
> 400 fim, respectively). Compression forces investigated are 
5, 10, 15, 20, 30, 40 and 50 kN, yielding sample thicknesses 
from 3.5 mm down to 2.7 mm. The samples are similar to 
those in Ref. [23 J . for which mercury intrusion porosimetry 
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showed than predominate pore size is around 1-5 |J.m. Lactose 
and calcium phosphate samples, on the other hand, were man- 
ufactured from pure powders with a 5-10 /im mean particle 
size. Compression forces were 30, 40 and 55 kN for lactose 
(thickness range 2.7-2.6 mm), and 20, 30, 40 and 50 kN for 
calcium phosphate (thickness range 2.1-1.9 mm). 

Diffusion modeling 

In a system with cylindrical symmetry, the time evolution 
of the intensity transmitted through a slab at a radial distance 
6r from the incident light pulse, /tC^, 5r), is given by Eq. [T| 

iHiia. 

There /q is the incident intensity, c the speed of light in 
the medium, the absorption coefficient, and and 
the radial and longitudinal diffusion constants respectively. 
The quantities z+.m and m contain the dependence on slab 
thickness L. More specifically, z^^^n = 2m{L + 2z^) + zq 
and z_,„i = 2m(L + 2z^) — 2ze + zo- Apart from z+,o = zq, 
which gives the position of the isotropic light source that mod- 
els injected light, they correspond to positions of negative and 
positive virtual sources introduced to handle boundary effects 
[^6|. The parameter Zg is the position of the, so called, ex- 
trapolated boundary. It is important to keep in mind that both 
zq and Ze depend on the diffusion constant D^- The source 
is located one transport mean free path into the medium, i.e. 
zq — 2>Dz/c, and Zg depends on the effective refractive index 
(rig) in a more complex manner. In fact, rig is not known ex- 
actly for our complex porous systems, but as long as sample 



dimension significantly exceed the transport mean free path 
of light, the time evolution of the transmitted intensity de- 
pends only weakly on zq and Zg. We use the approximate 
value Tig — 1.5 throughout this work, resulting in a Zg of 
3.6279 X 2Dz/c- Note that this value of Zg is calculated using 
ordinary isotropic considerations 1261 . and is only a conve- 
nient approximation for anisotropic media. We have verified 
that a change in the assumed value of rig has negligible influ- 
ence of results presented. It should also be noted that the size 
of the fiber cores (600 /im), i.e. the fact that 5r is not really a 
single value for a single measurement, are taken into account 
by integration over source and detector areas. Finally, note 
that the isotropic model of dynamic light diffusion through a 
slab is reached simply by setting = Dr = D, where D 
then is the isotropic diffusion constant of the medium. Fitting 
of Eq. [T] to experimental data is made using only diffusion 
constant(s), absorption coefficient and amplitude coefficients 
as free fit parameters (since zq and Zg depend on D^, they are 
varied accordingly). When multiple TOE histograms are eval- 
uated simultaneously (modeling of 5r series as one compound 
dataset), each single dataset has its own free amplitude coef- 
ficient. Measuring diffuse intensity in an absolute manner is 
extremely difficult, and attempts to maintain relative intensity 
relation for different Sr would be highly inefficient due to the 
enormous difference in transmission. The absolute time scale 
{t = 0) is inferred from careful measurements of the instru- 
mental response function (IRE), and the shape of the IRE is 
taken into account during fitting by measuring the fit in terms 
of agreement between experimental data and the convolution 
between IRE and diffusion model 1,25] . 
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